Anomalous proximity effect in an inhomogeneous disordered superconductor 
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By combining very low temperature scanning tunneling microscopy and spectroscopy on a TiN 
film we have observed a non uniform state comprising of superconducting (S) and normal (N) areas. 
The local density of states displays a spatial dependence between S and N different from the usual 
proximity effect. We conclude that mesoscopic fluctuations might play a major role in accordance 
with recent theories on superconductor-normal metal quantum transition. 

PACS numbers: 74.45.+C, 74.81.-g, 71.30.+h 
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In two dimensions, the superconductor-insulator tran- 
sition (SIT) is traditionally described by two microscopic 
mechanisms Q . In the first one , the Cooper pairs are all 
formed at T=0 and become localized through the SIT. 
This is often refered to as the "bosonic" model and has 
been successfully to granular superconductors Q The 
onset of the superconductivity appears at a constant tem- 
perature but the transition becomes broader as the sheet 
resistance Ra approaches jpj . This broadening is associ- 
ated with quantum fluctuations of the phase of the order 
parameter due to the competition between the charging 
energy of the superconducting grains and the Joseph- 
son coupling between them. The second mechanism at- 
tributes the weakening of the superconductivity to a dis- 
order enhanced Coulomb repulsion. Pair breaking is con- 
sidered simultaneously with the decreasing of both the 
superconducting critical temperature and the amplitude 
of the order parameter as the SIT is approached. This 
so-called "fermionic" scenario usually describes homoge- 
neous thin films which keep a sharp superconducting 
transition until disorder becomes very close to a critical 
value |3(. By further increasing the disorder, by reducing 
the thickness or by application of a magnetic field, these 
films can be driven into an insulating state with local- 
ized electronic excitations. However, real samples may 
be neither purely bosonic or purely fermionic systems |J| 
and a unified theory describing these two limiting cases is 
still needed. Moreover, there are theoretical predictions 
H and experimental observations 0, S of an inter- 
mediate metallic phase right at the SIT. 

In three dimensions, depending on the disorder, both 
a SIT and a superconductor-normal metal transition 
(SMT) are predicted at T=0 Q, 0] and have been exper- 
imentally observed in a field-tuned transition [lflj . The 
existence of an intermediate metallic phase is often ac- 
companied by a re-entrant temperature and field depen- 
dence of the electrical resistance R of the film . Ac- 
tually, the competition between Cooper pairs and quasi- 
particule tunneling between superconducting grains or 
clusters can explain a nonmonotonic behavior of R. In 
these models, a zero resitance state appears when the 
network of superconducting islands percolate [l^ . In- 
homogeneities in the morphology of granular films [l3| 



or strong statistical disorder in nominally homogeneous 
films can lead to the formation of these islands 
Mesoscopic effects induced by local fluctuations in the 
density of randomly distributed impurities can also drive 
the film into a multiple reentrant field-tuned SMT [l5| . 

A double reentrant behavior very close to a SIT has 
recently been observed in disordered TiN thin films . 
This SIT could be triggered either by a magnetic field or 
by increasing R n ■ It is worth noticing that such a double 
reentrance has also been observed in Josephson junction 
arrays and highly granular superconductors 0. In this 
Letter we report spectroscopic measurements on a simi- 
lar TiN film with a lower R a and in zero field, thus far 
away from the above mentioned SIT. Nevertheless, we 
show that disorder induced inhomogeneities already ex- 
ist in such a film which cannot be detected by transport 
measurements. 

TiN was prepared by DC reactive magnetron sputter- 
ing at 350°C on thermally oxidized Si substrates. By 
sputtering a Ti target at various nitrogen partial pres- 
sure in an argon-nitrogen gas mixture, different T1N5 
compounds with 0.7 < 5 < 1.2 can be obtained 0. We 
have been able to vary the room temperature electrical 
resistivity between 80 /ifi.cm and 1100 /iil.cm by chang- 
ing the nitrogen flow rate from 40 to 200 seem 0|. In 
stoichiometric TiN films, the mean free path I is mainly 
limited by the grain boundaries and one gets I ~ L g , 
where L g is the typical grain size. When the nitrogen 
flow rate is increased during the film preparation, Ti va- 
cancies are introduced and one obtains overstoichiomct- 
ric TiN<5 (8 > 1) with I < L g . Moreover the Ti vacan- 
cies are not uniformly distributed inside the grai ns but 
are rather concentrated at their boundaries 20]. The 
film studied here has an intermediate resisitivity of 270 
/ifLcm, a mean free path I = 5 ± 1 nm and a thickness 
of 100 nm 19]. We estimate L g ~ 20 nm from STM 
images. It undergoes a sharp superconducting transition 
at T c = 4.68 K as detected by R(T) experiments. 

Wc have combined topography and spectroscopy mea- 
surements with a scanning tunneling microscope (STM) 
cooled down to very low temperature in a dilution re- 
frigerator. In order to probe the local density of states 
(LDOS), a small AC modulation of 20 fiV rms is added 
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FIG. 1: Several spectra taken at different temperatures at 
the same location and the corresponding BCS fits. Inset : 
temperature dependence of the BCS gap 

to the sample-tip DC bias voltage V and the differen- 
tial conductance ^ is obtained with a lock-in amplifier 
technique. At any position above the sample V is ramped 
and the resulting ^ (V) curve gives the LDOS around 
the Fermi energy. Several spectra are displayed in FigQ] 
at different temperatures . The superconducting gap val- 
ues A(T) are obtained by convoluting a BCS density of 
states with a Fermi distribution function and are dis- 
played in the inset. A weakly temperature dependent 
Dynes parameter in the range T = 0.022mV to 0.027 mV 
is used to adjust the peak height j2l|. The lowest indi- 
cated temperature of 250 mK is not the measured one 
which was 125 mK but rather the temperature needed to 
correctly fit the spectrum. This indicates that the energy 
resolution of our STM is probably limited by unfiltered 
electromagnetic radiations which heat the electrons. We 
find A(T=0) = 0.73 mV and a ratio -£r = 1.81 not far 
from the theoretical BCS value of 1.76. 

However, these spectra are not observed everywhere 
on the surface of the TiN film and at other locations a 
normal-metal like flat LDOS is measured. In order to get 
an image of the superconducting and normal areas we set 
the bias to a voltage Vo = 0.75 mV slightly above the 
BCS gap value. Two images were then recorded simul- 
taneously : the topographic one Z(x,y) and the spectro- 
scopic one ^ (x,y) at a given energy eV . When the tip 
is scanning above a superconducting region, the differen- 
tial conductance signal increases because of the peak in 
the LDOS at Vo- Inversely, a flat LDOS is characterized 
by a lower output of the lock-in amplifier. The two im- 
ages can be merged into a 3D coloured picture as shown 
in FigEl where the dark orange areas are normal while 
the light yellow ones are superconducting 

Sometimes, isolated grains stand higher on top of the 
surface. They can be either normal or superconducting 
as is shown, respectively, in the upper and lower picture 
of Fig[3 Once the electronic properties are identified 




FIG. 2: 3D images of TiN : 400 x 400 nm 2 taken at 143 mK 
(top) and 250 x 250 nm 2 taken at 258 mK (bottom); the 
maximum height amplitude is 2nm for both pictures. The 
colouring maps the LDOS; the superconducting regions are 
yellow and the normal ones are orange. 



everywhere on the scanned surface, it is then possible to 
measure the complete spectra along a line which crosses 
the boundary between a superconducting and a normal 
region. Several spectra along such a scan are displayed in 
FigEH Except for the trully gaped one which is labeled as 
the starting position x exp — on the S side, none of the 
spectra are BCS-like. The distance dsN over which the 
LDOS is varying depends on the local granularity. It can 
be as long as 50 nm when the surface is smooth. But if 
we scan across an isolated superconducting grain such as 
those in the bottom panel of FigEl the spatial variation 
of the LDOS is much more rapid and takes place within 
less than 10 nm, mainly inside the grain itself. However, 
these different spatial scales do not modify the form of 
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FIG. 3: Spatial evolution of the LDOS between a super- 
conducting and a normal region. The dotted lines are ex- 
perimental data. The overall measured proximity length is 
dsN — 48 nm. The solid lines are numerical fits with Eq.© 
which has no adjustable parameter and the boundary condi- 
tions Eq.JJ). 



the overall transition of the LDOS at the SN interface 
since the spectra evolve in the same manner between S 
and N wherever they are obtained on the surface of the 
sample. 

The spatial dependence of the proximity effect at an 
SN interface can be described in the framework of the 
quasiclassical Green's functions by a complex pairing 
angle 9{E,r) The LDOS n at a position r and 

for an energy E from the Fermi level is related to 6 
by n(E,r) = noRe[cos9(E,r)]. In the dirty limit, i.e. 
I < £ = ^hD/2A (D is the diffusion constant), 6(E,r) 
obeys the Usadel equation which in one dimension can 
be written as|2J|: 



HD d 2 9 



[iE-T in - 2Y sf cos<9] sin6> 

+A(a;)cos6> = 0, 



(1) 



where T s f and Ti n are the spin-flip and the inelastic 
scattering rates respectively. This description of the su- 
perconducting proximity effect at a mesoscopic scale has 
already been checked experimentally in SN heterojunc- 
tions with nanofabricated tunnel junctions and with STM 
experiments [2^. For infinite SN systems, it has been 
shown that the LDOS exhibits a V-shaped pseudo-gap 
with peaks below and above the Fermi level. These peaks 
are separated by an energy which decreases as a function 
of the distance from the interface in the normal side of 
the junction. Here, we observe instead U-shaped spectra 
with peaks that are pinned at the BCS energy A for any 
position between S and N. The LDOS is flat for energies 
smaller than A and increases as the tip is moved progres- 
sively away from the S region. We tried to fit these results 
using Eq.JQ) and a self-consistently determined space de- 
pendent order parameter A(x) [26J. However, no set of 



parameters was able to correctly reproduce even qualita- 
tively the shape of the spectra. One of the possible rea- 
sons for this failure could be that our film does not fulfill 
the dirty limit condition. Actually, we have £ = 6.5 nm 
not much bigger than I 19]. Moreover, since the scat- 
tering centers are concentrated at the grain boundaries, 
the electronic trajectories could be considered to be qua- 
siballistic inside the grains. The projection in only one 
dimension of the Usadel equation could also be a too 
crude approximation. 

It is nevertheless striking than we can reproduce with a 
good accuracy our results if we assume the pairing angle 
varies linearly with distance : 
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integrated with the boundary conditions 
6(x = 0) = Obcs an .d 6(x 



and 9n = 0. Eq.(J2J has been 
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The numerical results are shown in FigOU We found an 
excellent agreement between the experimental data and 
the calculations without any adjustable parameters. A 
and r are deduced from the spectrum taken above the 
BCS superconducting region and positions are chosen in 
order to match the measured LDOS at the Fermi energy. 
The slight differences between these positions and the 
measured ones reflect the granularity always present; the 
latter affects the measured value of the proximity length 
ds 7v and can also perturb the pure linear behavior of 
Eq.(J2|. Although, there is no obvious physical meaning 
behind this equation, we want to point out the univer- 
sal character of this observed anomalous proximity effect 
since it does not depend on any physical properties of the 
material such as inelastic scattering times, conductivities 
or coherence lengths of S and N necessary to describe the 
usual proximity effect [13] • 

We want to discuss now the physical origin of the 
superconducting and normal clusters. In STM experi- 
ments, surface contamination must be considered seri- 
ously as a possible artefact for samples exposed to air. 
Nevertheless, there are easily recognizable signs which 
indicate the presence of adsorbates : (i) they usually de- 
grade the quality of the images, (ii) At moderately low 
temperature, when they are frozen out, their positions 
are revealed as bumps in the topographic images, (iii) At 
very low temperature, they generally exhibit the LDOS of 
an insulator and give dark areas in spectroscopic images. 
As seen in the pictures of FigO no such effects are visi- 
ble. The transition between superconducting and normal 
areas is in places very smooth with no change in the noise 
in the data. This proves that the inhomogeneities are in- 
trinsic to the sample. However, we cannot rule out that 
because of oxidation of intergrain boudaries, for exam- 
ple, the surface of TiN itself is not representative of the 



4 



bulk which could be more homogeneous. Another pos- 
sible source of inhomogeneities, could come from spatial 
variations of the nitrogen concentration in TiN^, since 
the BCS coupling constant, A depends on 5 [28j. Follow- 
ing Ioffe and Larkin's pionnering work, strong statistical 
fluctuations of A(r) can lead to an inhomogeneous system 
where superconductivity appears first in localized drops 
which percolate at T c even far away from the SIT thresh- 
old [14j. A spatial distribution of the order parameter 
amplitude would give a broadened density of states as 
obtained by Hsu et al with large area tunnel junctions 
Here, we would expect instead to observe well identi- 
fied superconducting clusters with different BCS gap val- 
ues. It seems therefore more likely that our film is made 
of well coupled superconducting grains with a unique A 
embedded in a normal metal matrix where the Ti vacan- 
cies pile up. According to Feigel'man et al and Spivak et 
al, quantum fluctuations in such a system make super- 
conductivity unstable [2fl| | . Spatial variations of the size 
and the concentration of grains as well as mesoscopic fluc- 
tuations of the intergrain conductance can drive the film 
locally into either a normal metal or a superconducting 
cluster. This is also consistent with the observed leveling 
off of the resistivity at very low temperature and with the 
reentrant field tuned transition observed in thinner TiN 
films. Indeed, when the film is closer to the percolation 
threshold of the superconducting network its global re- 
sistivity is then governed by very few bottlenecks and the 
associated mesoscopic fluctuations. According to Spivak 
et al, these mesoscopic fluctuations can give a multiple 
reentrant transition between superconducting and nor- 
mal metallic states |l5| . 

In summary, we have observed inhomogeneities in a 
disordered superconducting thin film of a nature simi- 
lar to those invoked to explain macroscopic properties 
000 and possibly described by recent fluctuations 
driven SMT models 29] . It is very likely that in order to 
understand the spatial dependence of the LDOS between 
superconducting and normal areas, one should go beyond 
the mean field Usadcl equation : mesoscopic fluctuations 
could help explain the inhomogeneous superconducting 
state, the re-entrant field tuned transition in TiN and 
modify the proximity effect 30]. 

We are grateful to M. Houzet for helping us in resolv- 
ing self-consistently Usadel equation and to M. Sanquer 
who initiated these experiments. 
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